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ABSTRACT 

The  generation  of  long-wavelength,  viscous-inviscid  interactive  Gortler  vortices  is  studied  in 
the  linear  regime  by  numerically  solving  the  time-dependent  governing  equations.  It  is  found 
that  time  dependent  surface  deformations,  which  assume  a  fixed  nonzero  shape  at  large  times, 
generate  steady  Gortler  vortices  that  amplify  in  the  downstream  direction.  Thus,  the  Gortler 
instability  in  this  regime  is  shown  to  be  convective  in  nature,  contrary  to  the  earlier  findings 
of  Kuban  and  Savenkov.  The  disturbance  pattern  created  by  steady  and  streamwise-elongated 
surface  obstacles  on  a  concave  surface  is  examined  in  detail,  and  also  contrasted  with  the  flow 
pattern  due  to  roughness  elements  with  aspect  ratio  of  order  unity  on  flat  surfaces.  Finally,  the 
applicability  of  the  Briggs- Bers  criterion  to  unstable  physical  systems  of  this  type  is  questioned 
by  providing  a  counterexample  in  the  form  of  the  inviscid  limit  of  interactive  Gortler  vortices. 
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1  Introduction 


The  physical  realizability  of  an  unstable  equilibrium  solution  corresponding  to  a  nearly  homo¬ 
geneous  shear  flow  is  dependent  on  whether  or  not  the  equilibrium  state  corresponds  to  the 
time-asymptotic  solution  of  the  related  unsteady  problem.  The  answer  to  this  latter  question 
can  usually  be  obtained  quite  easily,  provided  one  has  sufficient  information  to  classify  the  lin¬ 
ear  instabilities  of  this  equilibrium  state  as  being  of  either  absolute  or  convective  nature  (Briggs 
19(71,  Bers  1975,  Huerre  and  Monkewitz,  1985).  The  presence  of  an  absolute  instability  means 
that  any  compact  unsteady  source  with  a  continuous  Fourier  spectrum  will  introduce  distur¬ 
bances  that  amplify  indefinitely  in  time  at  all  neighbouring  locations,  thereby  eliminating  the 
possibility  of  any  time-asymptotic  state  within  the  linear  framework.  In  contrast,  the  response  of 
a  convectivelv  unstable  flow  to  a  pulsed  source  takes  the  form  of  an  instability-wave  packet  that 
is  swept  away  from  the  source  by  the  mean  flow,  thus  restoring  the  original  state  at  sufficiently 
large  times. 

For  convectively  unstable  flows,  it  is  also  possible  to  investigate  the  response  of  the  flow  to 
some  continuous  forcing  in  a  local  region,  say  in  the  form  of  a  time-harmonic  disturbance  which 
is  switched  on  at  a  finite  time.  The  time-asymptotic  solution  to  a  “signalling”  problem  of  this 
type  can  be  obtained  either  from  a  superposition  of  the  wave-packet  solutions  Balsa  (1988),  or 
by  directly  solving  the  steady  state  equations  which  have  the  harmonic  time  dependence  built 
into  them.  The  overall  amount  of  effort  required  is  perhaps  larger  in  the  former  case;  however, 
it  automatically  ensures  the  causality  of  the  steady-state  response,  without  any  need  (as  in  the 
latter  rase)  for  extraneous  considerations  such  as  the  Briggs-Bers  criterion  in  order  to  specify 
the  correct  streamwise  boundary  conditions  at  the  upstream  and/or  downstream  ends. 

The  objective  in  this  paper  is  to  consider  the  above  issues  in  the  context  of  a  centrifugal,  i.e., 
the  Gortler  instability  of  boundary-layer  flows  over  concave  surfaces.  Experiments  have  shown 
that  this  instability  is  usually  manifested  in  the  form  of  steady  streamwise  vortices  which  amplify 
in  the  downstream  direction,  eventually  yielding  to  time-dependent  instabilities  of  a  secondary 
and  tertiary  nature.  The  various  experimentally  observed  features  of  the  linear  and  nonlinear 
development  of  the  Gortler  vortices  have  also  been  explained  using  asymptotic  theories  (Hall 
1982a.-b,  Hall  and  Lakin  1988,  Hall  1990),  and  direct  numerical  simulations  (Hall  1988,  Sabry 
and  Liu  1990,  Liu  and  Sabry  1991). 

While  the  occurence  of  steady  unstable  vortices  is  suggestive  of  a  convective  nature  for  the 
Gortler  instability,  the  validity  of  this  assumption  had  not  been  addressed  until  Parks  and  Huerre 
(1988),  who  examined  the  special  case  of  an  asymptotic  suction  profile.  More  recently,  Ruban 
(1990)  (herafter  referred  to  as  R)  examined  the  Gortler  instability  of  an  arbitrary  profile  in  the 
long-wavelength  regime  using  asymptotic  methods  in  the  limit  of  a  large  Reynolds  number  /?,  and 
a  large  Gortler  number  G.  Somewhat  surprisingly,  he  found  that  the  flow  is  absolutely  unstable 
to  spanwise-periodic  perturbations.  However,  he  also  examined  the  disturbances  created  by  an 
isolated  point  source  and  concluded  that  the  flow  is  convectively  unstable  in  this  latter  case. 
Here,  the  wave  packet,  determined  through  steepest-descent,  analysis,  was  found  to  spread  in  all 
directions  including  upstream,  except  for  a  sector  of  60°  directly  ahead  of  the  source. 

Savenkov  (1990)  (henceforth,  S)  applied  the  same  asymptotic  framework  as  Ruban  to  study 
the  receptivity  problem  for  the  Gortler  vortices  in  the  long- wavelength  regime.  He  used  a  residue 
analysis  in  conjunction  with  Fourier  transform  methods  to  compute  the  approximate  temporal 
development  of  the  perturbations  created  by  an  unsteady  wall  hump  which  assumes  its  final 


(nonzero)  shape  after  a  finite  length  of  time.  His  results  were  supportive  of  the  conclusions 
in  R.  in  that  the  steady-state  pattern  was  found  to  he  similar  in  shape  to  the  time-dependent 
wave  packet  of  R.  and  no  (Idrtler  vortices  were  generated  downstream  of  the  hump.  Previously. 
Rozhko  (t  al  (  HISS)  (henceforth,  RRT)  had  also  studied  the  same  problem,  albeit  by  solving  the 
steady  state  equations  directly  in  Fourier  transform  spare.  The  disturbance  pattern  computed 
in  this  manner  was  also  similar  to  that  obtained  by  S  as  a  limit  of  the  time-dependent  problem. 
RIFF  had  also  found  that  a  crucial  difference  between  the  disturbance  patterns  on  concave  and 
convex  walls  corresponds  to  the  presence  of  an  0(1)  upstream  signature  in  the  former  cast',  as 
against  a  total  lack  of  any  upstream  influence  in  the  latter  case. 

The  above  results  concerning  the  absolute  nature  of  the  Oort ler  instability,  as  well  as  the 
occurence  of  an  upstream  interaction  in  the  concave-wall  case,  are  rather  intriguing,  the  former, 
since  it  casts  a  doubt  on  the  relevance  of  much  of  the  previous  work  on  (idrtler  instability  which 
was  related  to  steady  vortices,  and  the  latter,  due  to  the  fact  that  the  steady-state  equations 
are  actually  parabolic  in  nature.  In  this  paper  we  examine  the  problems  addressed  by  R,  S  and 
RRT,  and  demonstrate  that  the  peculiar  findings  of  these  investigators  are  actually  a  result  of 
having  ignored  a  fundamental  property  of  the  part ial  differential  equations  governing  the  steady 
and  unsteady  problems.  Specifically,  we  show  that  solving  the  unsteady  problem  with  the 
appropriate  boundary  conditions  not  only  confirms  the  convective  nature  of  (Idrtler  instability 
in  the  spanwise-periodic  case,  but  also  yields  the  solution  to  the  receptivity  problem  in  the  limit 
of  large  t imes. 

The  problem  of  the  generation  and  subsequent  linear  development  of  long-wavelength  Gbrtler 
vortices  is  formulated  in  the  following  section.  Although  no  extra  effort  is  involved  in  solving 
the  complete  problem  related  to  viscous  inviscid  interactive  vortices,  it  is  more  illustrative  to 
separately  consider  the  two  limiting  case's  corresponding  to  the  viscous  and  inviscid  vortices, 
respectively.  Both  analytical  and  numerical  methods  an'  used  towards  this  purpose;  the  analyt¬ 
ical  work  is  presented  in  Section  2.  whilst  in  Section  '■]  we  present  our  results.  Finally  we  draw 
some  conclusions  in  Section  1. 


2  The  governing  equations  for  viscous-inviscid  interac¬ 
tive  longitudinal  vortex  structures  and  some  analyti¬ 
cal  results 

At  a  given  value  of  the  curvature  parameter,  riz.  the  (Idrtler  number  (7  =  ^7 /?c^.2,  the  linear 
development  of  steady  (Idrtler  vortices  is  determined  primarily  by  the  nondimensional  wave¬ 
length  in  the  spanwise  direction.  A-  =  li  c J/.2  jz  ■  Here  /,*  denotes  a  typical  streamwise  length 
scale  based  on  the  distance  between  the  leading  edge  and  the  location  of  interest,  n‘  (>>  0(L*)) 
is  the  radius  of  curvature  of  the  surface  at  this  location,  and  Rr /,•(>>  1 )  is  the  Reynolds  number 
based  on  the  local  free-stream  speed  (  *  ,  and  the  distance  /,*.  The  asterisk  is  used  to  indicate 
the  dimensional  quantities  here,  and  throughout  the  rest  of  ibis  paper. 

It  was  shown  independently  by  Denier.  Hall  and  Seddougui  (1990),  and  I  imoshin  ( 1990) 
that  in  the  limit  of  large  (Idrtler  numbers  ((»’  >>  1)  one  needs  to  consider  five  different  asymp¬ 
totic  regimes  along  the  A,  axis.  In  the  order  of  increasing  spanwise  wavelength,  these  regimes 
correspond  to  (i)  the  neutral  regime  (A-  ~  (/”''  *)■  (ii)  the  most  unstable  regime  (A-  ~  (1  1/S)- 


(iii)  the  inviscid  regime  (A-  =  0(1)),  (iv)  the  long  wavelength  or  viscons-inviscid  interactive 
regime  (Ac  ~  (71/').  and  finally,  (v)  the  nonparallel  regime  corresponding  to  A.  =  ()((i]^2). 
Investigations  of  the  linear  and  nonlinear  stability  problems  associated  with  the  first  three  of 
these  five  regimes  have  been  given  by  Hall  (I982a-b),  Hall  and  Lakin  (1988),  Denier,  Hall  and 
Seddongni  (1990),  and  Timoshin  (1990).  The  viscons-inviscid  regime  for  (Jortler  vortices  was 
first  investigated  by  Rozhko  and  Rnban  (1987)  in  connection  with  the  disturbances  generated 
bv  streamwise-elongated  roughness  elements  over  curved  surfaces,  and  subsequently,  in  the  same 
context  by  RRT,  R  and  S,  as  mentioned  already.  The  fifth  regime  corresponding  to  nonparallel 
vortices  does  not  appear  to  have  been  studied  in  any  detail  until  now,  but  see  H all ( 1 983 ) . 

As  stated  in  the  Introduction,  our  concern  in  this  paper  is  with  the  generation  and  linear 
amplification  of  the  long  wavelength,  or  viscons-inviscid  interactive  (Jortler  vortices,  which  are 
relevant  at  locations  close  to  the  leading  edge.  As  first  shown  by  RR,  the  asymptotic  scalings 
of  this  regime  are  fixed  by  the  condition  of  viscons-inviscid  interaction,  plus  a  balance  of  tlx* 
displacement-induced  pressure  with  the  jump  in  pressure  across  the  main  part  of  the  boundary 
layer  due  to  centrifugal  effects  associated  with  the  surface  curvature.  These  balances  identify 
the  spanwise  length  scale  A-.  based  on  the  boundary-layer  thickness,  as  being  of  0(<~').  when' 
<  —  (7-|A  (<<  1),  whereas  the  streamwise  length  scale  (s,  based  on  the  distance  L‘  from  the 
leading  edge,  is  determined  to  be  (7(f3).  In  the  normal  direction,  the  vortex  structure  consists 
of  three  distinct  subregions,  viz.  the  lower,  middle  and  upper  decks,  whose  thicknesses  relative' 
to  that  of  the  boundary  layer  are  of  0(f),  0(1),  and  0(f  ~ 1 )  respectively.  The  overall  dynamics 
of  each  of  the  decks,  as  well  as  their  coupling,  is  quite  analogous  to  that  in  the  conventional 
three-dimensional  triple  deck  problem.  The  structural  similarities  between  the  two  problems 
lead  one  to  anticipate  that  the  (Jortler  vortex  problem  is  also  amenable  to  the  same  solution 
procedure  as  that  applied  in  triple-deck  problems,  and  indeed,  this  was  shown  to  be  the  rase  by 
RR.  Specifically,  it  is  possible  to  obtain  closed  form  solutions  for  the  perturbations  in  the  middle 
and  upper  decks,  and  a  matching  of  these  two  provides  the  interactive  relationship  between  the 
pressure  and  displacement-thickness  perturbation,  this  closes  the  lower-deck  problem  governed 
by  some  form  of  the  three-dimensional  boundary  layer  equations. 

for  the  sake  of  definiteness,  consider  the  same  problem  as  that  studied  In  both  R  and  S. 
viz ..  the  disturbance  pattern  produced  by  an  unsteady  wall-hump  with  an  arbitrary  but  specified 
shape  that  can  be  expressed  as  the  coordinate  surface 


Y  =  0, 


(2.1) 


corresponding  to  the  Prandtl-transposed  coordinate  within  the  lower  deck,  V  -  ,  1  Hr) y  A4/1 1/*/ L~  - 
h  F(.Y,  Z,  7  ).  where  .Y  =  (-3A and  Z  =  (  He J(2A10F -*//,’  are  local  coordinates  in  the 
streamwise  and  spanwise  directions  respectively,  T  —  (2 //cj(2 A~*( '  I  ,*J* / L‘  is  the  nondimen- 
sional  time,  and  h  (<<  1)  and  F(A,  Z,  7  )  are,  respectively,  the  normalized  height  parameter 
and  shape  function  characterizing  the  wall-hump  geometry  (see  Fir.  1).  For  reasons  which  will 
be  obvious  later,  it  is  appropriate  to  only  consider  the  type  of  l  ump  geometries  which  have  a 
definite  streamwise  origin,  i.e.,  F(A  ,  Z,  7  )  =  0  for  A  less  than  some  finite  value  Ao>  which  will 
be  assumed  to  be  the  origin.  A  =  0.  without  any  loss  of  generality.  The  parameter  A  is  used 
to  denote  the  wall  shear  corresponding  to  the  incoming  boundary-layer  profile  just  upstream  of 
t  he  roughness  element . 


Introducing  the  expansions 


^:-)  =  (Y  +  f/,A-’/7F. 


Rf 


-1/2 


h \-'/7(Fr  + 


V  Fx  +  V),  f  ReLV2hX1! ‘  IT,  e~6 Ftcjl  h.\4/ '  P)  for  the  lower  deck  variables,  one  finds  that 
the  leading-order  perturbations  ((/,  V ,  W,  P)  are  governed  by  the  linearized,  three-dimensional 
boundary-layer  equations 

d\U  dy  V  +  dzW  =  0,  (2. 2a) 

(dr  +  Ydx)U+  V  =  d$lJ,  (2.2  b) 

and 

(dr  +  V'  dx)\V  =  -dzP  +  dyW,  (2.2c) 

without  any  pressure  gradient  along  the  streamwise  direction.  Since  the  surface  deformation 
due  to  the  unsteady  hump  has  been  assumed  to  originate  at  a  finite  time,  the  disturbance  field 
may  be  taken  to  be  zero  at  the  initial  instant  of  time.  Similarly,  it  can  be  assumed  that  the 
surface  obstacle  will  not  produce  any  disturbances  sufficiently  far  upstream,  and  also  that  the 
disturbance  motion  either  decays,  or  remains  bounded,  as  \Z\  — >  oo,  depending  on  whether 
F(.X ,  Z,T)  has  a  compact  or  noncompact  support  in  the  spanwise  direction.  The  boundary 
conditions  in  the  normal  direction  are  given  by 

V  =  V  =  W  =  0  at  Y  =  0,  (2.3a  -  c) 

and 

V  ->  F{ X,  Z,  T)  —  A(X,  Z,  T),  W  -♦  0  as  Y  ->  oc,  (2.3d,  e) 

wliere  the  displacement-thickness  perturbation  A( X,  Z,  T)  is  related  to  the  pressure  perturbation 
P(.X ,  Z.T)  via  the  interactive  relationship  (RR) 


SzP  =  -.-is m(G)c)zA  +  -  I*™  d*AiX7i]d(,  (2.3/) 

which  reflects  the  simultaneous  balance  between  the  lower-deck  pressure,  the  curvature-induced 
pressure  within  the  middle  deck  and  the  displacement-induced  pressure  inside  the  upper  deck. 
In  ( 2 . 3 f ) .  the  sign  of  the  Gbrtler  number  has  been  assumed  to  be  positive  for  concave  surfaces, 
and  negative  for  convex  ones.  Moreover,  one  may  also  obtain  results  for  the  case  of  an  underlying 
surface  that  is  flat,  or  has  an  asymptotically  small  curvature,  by  setting  sign(G)  =  0  in  (2.3f). 
In  that  rase,  the  small  parameter  (  is  to  be  interpreted  as  a  measure  of  the  aspect  ratio  of  the 
plauform  of  the  surface  obstacle. 

It  is  worth  mentioning  that  apart  from  the  absence  of  a  streamwise  pressure  gradient  term  in 
equations  (2.2a-c).  which  follows  as  a  consequence  of  the  length-scale  disparity  in  the  streamwise 
and  spanwise  directions,  the  problem  of  viscous-inviscid  interactive  vortices  differs  from  the 
three-dimensional  triple-deck  problem  in  two  other  aspects  which  are  related  to  properties  of  the 
middle  and  tipper  decks,  and  thus,  are  manifested  through  the  interactive  relationship  (2.3f). 
Firstly,  the  role  of  the  main  deck  is  no  longer  just  a  passive  one,  corresponding  to  a  direct 
transmission  of  the  outer  pressure  to  the  lower  deck,  since  the  centrifugal  effects  in  the  middle 
deck  now  substantially  alter  the  pressure  gradient  imposed  on  the  lower  deck.  As  seen  from 
equation  (2.3f).  the  pressure  jump  across  the  middle  deck  turns  out  to  be  independent  of  the 
details  of  the  incoming  boundary-layer  profile,  being  equal  to  just  the  normalized  displacement- 
thickness  A  in  terms  of  its  magnitude.  For  this  reason,  the  equation  set  (2. 2-2. 3)  is  a  canonical 
one  and,  indeed,  after  some  additional  analysis,  it  can  be  shown  to  be  valid  at  both  subsonic  and 
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supersonic  speeds.  The  other  source  of  differences  between  the  Gortler-vortex  and  triple-deck 
problems  is  related  to  the  behaviour  in  the  upper  deck  where,  due  to  the  negligible  pressure 
gradient  in  the  streamvvise  direction,  the  secondary  flow  in  the  cross-flow  plane  is  decoupled 
from  the  local  streamvvise  motion,  in  a  manner  somewhat  analogous  to  slender  body  theory.  As 
a  result,  the  upper  deck  motion  is  elliptic  only  within  the  cross-flow  (V-  —  Z )  plane  and,  thus, 
cannot  exert  an  0(  1)  upstream  influence  on  the  lower  deck  motion  via  the  pressure-displacement 
relation  (2.3f). 

Applying  a  Fourier  transform  in  the  spanwise  direction  ( Z  — >  A-),  and  denoting  the  trans¬ 
formed  variables  using  an  overbar,  the  governing  equations  can  be  rewritten  as 


dxU  +  dyV  +  ikW  =  0, 

(2.4a) 

(<h  +  Ydx)Cr  +  V  =  dyO, 

(2.4/,) 

and 

(dT  +  YdxW  =  -ilcP  +  4  ILL 

(2.4c) 

along 

with  the  boundary  conditions 

II 

*7' 

II 

II 

o 

II 

o 

(2.5  a  —  r) 

V  =  F( A\  Z,  T)  -  A(X,  Z,  T).  W  =  0  as  Y  -*•  oo. 

(2.5 d,  e ) 

and 

kP  =  —sign(G)kA  +  sign(k)  Axx- 

(2.5/) 

obtained  from  (2.3a-f),  plus  the  homogeneous  initial  conditions  in  A'  and  T. 

In  the  following  section,  we  present  numerical  solutions  to  the  unsteady  problem  posed  above 
for  the  case  of  spanwise-periodic  perturbations,  and  confirm  the  lack  of  any  absolute  instability 
by  demonstrating  the  appearance  of  an  unstable,  spatially-growing  Gortler  vortex  at  large  times. 
It  is  also  possible  to  obtain  closed-form,  analytical  solutions  by  taking  the  appropriate  transform 
along  the  streamwise  direction.  For  instance,  RRT,  R  and  S  used  a  two-sided  Fourier  transform 
in  their  work;  but,  one  should  note  that  using  a  Fourier  transformation  along  the  streamwise 
direction  in  unstable  physical  systems  is  not  always  a  straightforward  matter,  although  fairly 
standardized  recipes,  such  as  the  Briggs-Bers  criterion  (see  Briggs  1964,  Bers  1975),  are  now 
available  to  treat  this  issue  in  the  case  of  streamwise  elliptic  systems.  Now,  due  to  a  lack  of  the 
streamwise  pressure  gradient  as  well  as  streamwise  diffusion  terms,  the  system  (2. 2-2. 3)  is  not 
elliptic  in  A  and  V’.  In  fact,  in  the  limit  of  very  long  wavelength  vortices,  it  can  be  easily  shown 
to  be  parabolic  in  the  streamwise  direction;  refer  to  the  discussion  following  Eq.  (2.1 2d).  Noting 
the  failure  of  RRT,  R,  and  S  in  accounting  for  this  anticipated  lack  of  upstream  influence  in 
the  concave-wall  case,  it  would  seem  that  a  straightforward  application  of  the  Fourier  transform 
technique  in  noil-elliptic  unstable  systems  may  lead  to  physically  unacceptable  results;  In  Section 
3.1  below,  we  will  also  consider  the  short-wavelength  (or  inviscid)  limit  of  the  problem,  and 
illustrate  the  inapplicability  of  Briggs-Bers  criterion  towards  predicting  the  nature  of  instability 
in  that  case. 

Thus,  consider  the  Laplace  transform  solution  to  the  problem  posed  by  equat  ions(2.4-2.5). 
which  is  given  by 


^Ai'(Co)  +  K«/»(<7)  A-2  -  |A-M2]  £  Ai(C)  d( ' 


(2.6fl) 


kP  =  [//2.s  ign(k)  —  ksign(G)]A, 


(2.6b) 


'  1W  +  116,(0  '^o)^011’ 


(2.6c) 


(2.6rf) 


and  V’  follows  from  the  A’-momentum  equation,  (2.4b).  Here  cr  and  /)  denote  the  transform 
variables  corresponding  to  the  time  T.  and  the  streamwise  coordinate  A  respectively,  while  (  is 
defined  as 


C  =  J  +  Co.  Co  = 


(2.6c,/) 


The  dispersion  relationship  corresponding  to  the  Cortler  vortex  instability  modes  is  then  given 
bv 


DW.  .1: k)  =  +  «XG)*2  -  \k\D2  =  0. 

MO  A 


(2.7 a) 


( L )  (A/)  (V) 

which  further  reduces  to 

:M>'(0)/5/3  +  sign(Ci)k2  -  \k\(i 2  =  0,  (2.7 b) 

for  the  case  of  steady  vortices  (a  =  0).  which  are  apparently  more  important  in  practice  than 
the  unsteady  ones.  An  interesting  property  of  the  steady  dispersion  relation  (2.7b)  is  that  it 
possesses  a  unique  root  in  the  complex  /  plane  for  any  given  (real)  spanwise  wavenumber  k\  and 
moreover,  this  root  always  lies  on  the  real,  positive  /  axis,  implying  that  the  steady  vortices 
exhibit  a  purely  exponential  growth  in  the  streamwise  direct  ion.  A  plot  of  the  stationary  growth 
rate  S  as  a  function  of  the  spanwise  wavenumber  k  is  shown  by  the  solid  curve  in  Fig.  2,  where 
the  large  and  small  k  asymptotes,  given  by  (2.8b)  and  (2.9b)  below',  are  also  indicated  by  the 
two  dashed  curves. 

As  described  earlier,  both  the  unsteady  as  well  as  steady  dispersion  relations  represent  a 
balance  between  the  effects  of  viscosity  in  the  lower  deck  (L),  destabilizing  centrifugal  forces 
in  the  middle  deck  (M),  and  the  viscous-inviscid  interaction  via  the  upper  deck  (U).  However, 
in  the  limit  of  large  spanwise  wavenumbers  (k  >>  1),  and  commensurably  larger  frequencies 
(|(t|  >>  1).  the  dominant  balance  shifts  to  just  the  inviscid  terms  (M-U)  in  (2.7a,b)  together 
with  the  large  cr  form  of  (L).  and  is  given  by 


—cr.i  +  si<ju(Ci')k2  —  |A'|/2  =  0, 

in  the  unsteady  case.  For  the  steady  problem  (L)  is  negligible  and  we  obtain 


(2.8a) 


s><l>)((i')k2  —  |Aj//2  =  0. 


(2.86) 
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which  matches  with  the  dispersion  relation  in  the  main  inviscid  regime  corresponding  to  sta¬ 
tionary  vortices  with  spanwise  wavelengths  comparable  to  the  boundary-layer  thickness;  see 
Denier,  Hall  and  Seddougui  (1990),  and  also  Timoshin  (1990).  In  contrast,  the  opposite  limit 
of  extra-long  spanwise  wavelengths  (k  <<  1),  which  was  the  focus  of  the  work  by  RRT,  R  and 
S,  leads  to  the  purely  noninteractive,  viscous-centrifugal  (L-M)  balances 


«„)  .  ,  „„,2  „ 

~  T  .S/-<77/((  »  )ft  —  0, 

f£Ai(C)dc 

(2-9  a) 

3/W'(0)/*5/3  +  .s  ign{G)k2  =  0, 

(2.9ft) 

in  the  unsteady  and  steady  cases,  respectively. 

Finally,  setting  sigv(Cl)  =  0  in  (2.7a)  yields  the  interactive,  flat-surface  limit, 


03Ai'«;o) 

J£Ai(Qdt 


-  W1  =  o, 


(2.10) 


which  corresponds  to  very  oblique  ToIImien-Schlichting  waves  (Hail  and  Morris  1991),  and  could 
also  have  been  obtained  by  taking  the  appropriate  limit  of  the  full  three-dimensional  triple  deck 
equations.  Since  our  primary  interest  lies  in  the  CJdrtler  vortex  type  of  instabilities,  the  limits 
(2.8)  and  (2.9)  are  more  relevant  to  us  than  (2.10),  and  therefore,  we  shall  choose  to  concentrate 
on  these  two  cases  in  the  remaining  part,  of  this  paper.  Of  particular  interest  will  be  the  extra- 
long  wavelength  problem  which  is  exemplified  by  the  dispersion  relationships  (2.9a,b)  and,  was 
also  considered  by  RRT,  R,  and  S.  For  this  case,  it  is  possible  to  further  substantiate  the  lack 
of  any  upstream  influence  in  a  manner  described  below. 

Basically,  further  manipulation  of  equations  (2.4)  leads  to  a  single  partial  differential  equation 
for  the  vertical  velocity  perturbation, 


\i%-(<h  +  YdxWV  =  0, 

(2.11) 

which,  toget  her  with  the  boundary  conditions 

V  =  ()y  l  =  0  at  V  =  0, 

(2.1 2a  -  b) 

f)yV(oo)  =  0, 

(2.12c) 

and 

(h()y  V'(0)  +  */>(< 7)  k2  [rVV'(oo)  +  i)x  F(  .V.  Z.  7’)]  =  0. , 

(2.12c/) 

and  the  homogeneous  initial  conditions  in  X  and  7  directions,  can  be  viewed  as  a  one  parameter 
(  viz.  )  )  family  of  linear  ad  vert,  ion  equations  in  the  ,Y  —  T  space  with  a  positive  semi-definite 
range  of  convection  velocities,  the  member  equations  being  coupled  via  normal  diffusion,  and  the 
integral  evolutionary  constraint  (2. Rid).  1  he  nature  of  these  advert. ion  equations  lends  further 
support  to  our  previous  argument  concerning  the  lack  of  any  upstream  influence  in  the  problem. 
Since  the  propagation  of  disturbances  in  the  upstream  direction  is  prohibited,  the  homogeneous 
initial  condition  in  the  streamwise  direction  can  be  imposed  just  upstream  of  the  origin  of  the 
surface  nonuniformity,  i.e.  at  A  — *  0  —  .  One  should  also  note  that  the  above  classification  of 


the  governing  partial  differential  equations  is  independent  of  the  sign,  or  the  magnitude,  of  the 
curvature  parameter  C.  since  it  only  appears  through  the  coefficient  of  a  zeroth-order  derivative 
term  in  (2.12(1).  The  only  difference  between  the  convex  and  concave  cases  corresponds  to  the 
existence  of  a  (Ibrtler  instability  in  the  latter  case,  which  leads  to  disturbance  amplification  in 
the  downstream  direction.  This  aspect  appears  to  have  been  overlooked  by  RRT,  R  and  S.  who 
utilized  a  two-sided  Fourier  transform  in  the  streamwise  direction  without  accounting  for  the 
properties  of  the  governing  differential  equations. 


3  Numerical  Results 

In  this  section  we  present  the  numerical  solutions  to  the  set  of  governing  equations  corresponding 
to  (2.2)  and  (2.2).  Tin*  cast'  of  spanwise  periodic  perturbations  is  examined  first  in  Section  3.1, 
followed  by  the  cast'  of  an  isolated  surface  excrescence,  which  is  discussed  in  Section  3.2.  In  both 
cast's,  spectral  discretizations.  Chebyshev  and  Fourier,  respectively,  were  utilized  along  the  V  and 
Z  directions,  whereas  a  second-order  accurate  backward  difference  scheme  was  used  to  compute 
tin'  vortex  evolution  in  space  (.V)  and  time  ( 7 ’ ) .  Since  both  the  l  and  IT  perturbations  approach 
their  limiting  values  at  infinity  (corresponding  to  equations  2. tube)  rather  slowly,  the  usage  of 
higher  order  boundary  conditions  obtained  through  asymptotic  considerations  was  found  to  be 
imperative  in  maintaining  the  spectral  accuracy  in  the  V’-direction.  Overall,  this  numerical 
scheme  combines  a  robust  marching  procedure  with  a  strong  coupling  in  the  secondary-flow 
plane,  the  latter  being  especially  desirable  for  solving  vortex-flow  problems.  The  accuracy  of 
the  numerical  results  was  verified  by  grid-resolution  checks,  as  well  as  through  comparisons  with 
analytical  solutions,  as  discussed  below. 

3.1  Unsteady  Evolution  of  Spanwise  Periodic  Perturbations 

We  will  first  discuss  the  results  pertaining  to  the  limit  of  tht'  extra-long  wavelength  viscous 
vortices.  We  recall  that  this  cast'  has  been  considered  by  R.  As  pointed  out  in  that  paper, 
and  as  can  be  seen  directly  from  tht'  dispersion  relation  (2.9a).  this  regime  admits  a  similarity 
behaviour  of  the  type  .V  ~  Z,i//r>  and.  T  ~  Z,/M.  anti  hence,  it  is  sufficient  to  consider  the 
unsteady  evolution  of  just  a  single  Fourier  mode  in  the  spanwise  direction;  this  we  take  to  be 
F  =  1  without  any  loss  in  generality.  The  real  and  imaginary  parts  of  the  spatial  (3)  root  of 
this  unsteady  dispersion  relation  art'  shown  in  Fig.  3  for  real  values  of  the  frequency  u.’(=  in). 
Fnlikt'  the  case  of  cross-flow  vortices,  where  the  most  unstable  modes  are  usually  time  dependent, 
here  the  steady  (Ibrtler  vortex  (a.*  =  0)  is  found  to  have  the  largest  growth  rate  for  any  given 
spanwise  wavenumber.  This,  along  with  the  fart  that  a  class  of  potentially  dominant  receptivity 
mechanisms  preferentially  excite  the  steady  vortices  as  compared  to  unsteady  ones  (Choudhari 
and  Streett,  1990).  would  help  explain  the  observed  dominance  of  steady  (Ibrtler  vortices  in 
laboratory  experiment s. 

As  shown  in  Fig.  3.  tht'  spatial  growth  rate  decreases  monotonically  with  frequency.  A 
neutral  point  exists  at  uv  =  2.29S  whilst  in  addition  there  is  the  asymptote  R<{3)  ~  as 

— >  x.  with  <l2 ( ~  — 3.2-1S2)  being  the  second  zero  of  A>'((0)  on  the  negative  real  Co  axis.  On 
the  other  hand,  the  imaginary  part  of  3  is  zero  al  ^  =  0.  implying  a  purely  exponential  growth 
by  the  stationary  vortex,  as  mentioned  previously;  but  Iw{3)  increases  nearly  monotonically  as 


S 


u?  is  increased,  leading  to  an  asymptotic  behaviour  of  lw(,i)  ~  J/^2  at  large  frequencies.  One 

may  also  observe  that,  since  the  growth  rate  curve  is  locally  stationary  near  its  maximum  at 
=  0,  the  derivative  c,,x  =  corresponding  to  th<’  group  velocity  of  the  stationary  vortices 
in  the  st Teamwise  direction,  is  purely  real  and  positive,  equal  to  about  1.5  as  seen  from  the 
dash-dot  curve  in  Fig.  3.  Of  course,  the  behaviour  of  Odrtler  vortices  is  far  from  that  of  the 
mainly-oscillatory  instability  modes  of  the  Tollmien-Schlicht ing  or  Rayleigh  type,  and  hence, 
the  notion  of  group  velocity  is  not  expected  to  be  physically  relevant  in  this  cast1,  especially  for 
the  u»'  =  0(1)  modes  when*  the  real  and  imaginary  parts  of  c,,\  are  comparable  to  each  other. 
Nevertheless,  it  is  interesting  to  note  that  the  real  part  of  the  group  velocity  in  the  X-direction 
is  always  positive,  and  less  than  its  value  for  stationary  vortices.  Furthermore,  the  imaginary 
part  of  the  group  velocity  is  bounded  for  all  spatial  modes.  Therefore,  one  might  expect  an 
impulsive  source  to  generate  a  wavepacket  with  both  ends  propagating  at  finite  speeds.  This 
was  indeed  found  to  be  the  case  in  the  numerical  solutions,  which  are  presented  below. 

For  t  he  purpose  of  simulating  t  he  t  ransient  wavepacket  problem,  we  considered  a  hump  shape 
of  t  he  form 


F(  .V,  Z,  7)  =  (1  -  cos  ( 2  7r  A  ) )  (1  —  cost  \kT))  cos(‘2kZ)  R(X.O. \)R{J\  0.0.5).  (3.1  a) 

when*  R(q.  </,.  r/f)  denotes  the  restriction  operator  in  the  q  space,  being  equal  to  unity  for  </,  < 
q  <  if/,  and  zero  elsewhere.  Similarly,  in  order  to  obtain  the  steady  hump  solution  in  a  causal 
manner,  an  unsteady  hump  shape  of  the  form 

F(.WZ.T)  =  (1  -  co.s(2tt.V))  (I  -  rr:'))  cos(2ttZ)  R(X. 0.1),  (3.1ft) 


was  also  considered. 

For  both  these  geometries,  we  have  shown  the  behaviour  of  the  displacement  function 
A(X,Z,  /  )  along  the  vortex  boundary  Z  =  0  as  a  function  of  .V  at  selected  instants  of  time  in 
Figs,  fa  and  lb  respectively.  Figure  (  la)  shows  that  the  transient  surface  deformation  generates 
a  vortex  patch  downstream  which  increases  in  its  streamwise  extent,  and  in  amplitude,  with  an 
increase  in  time.  One  may  note  that  between  0  <  A  <  1.  the  displacement  function  conforms 
quite  closely  to  the  hump-shape  during  the  period  0  <  7  <  1  when  it  is  nonzero.  Significantly 
the  vortex  pattern  displays  no  oscillations  in  the  streamwise  direction  at  all.  Thus,  the  bound¬ 
ary  layer  experiences  a  positive  displacement  at  all  times  along  the  vortex  boundary  Z  =  0. 
and  a  negative  displacement  along  the  centreline  Z  =  tc .  As  a  result  of  the  amplification  with 
the  passage  of  time,  one  would  observe  an  increasingly  stronger  streaky  structure  corresponding 
to  an.  alternate  pattern  of  boundary  layer  upliftment  (i.e.  deceleration)  and  an  equally  strong 
attachment  (i.e.  acceleration)  with  a  constant  spacing  along  the  streamwise  direction.  It  will, 
of  course,  be  very  interesting  to  see  how  nonlinearity  modifies  this  behaviour,  and  this  will  lx* 
the  topic  of  a  forthcoming  paper.  Finally,  but  most  importantly,  the  convective  nature  of  tin' 
spanwise  periodic  perturbation-;  is  seen  very  clearly  from  the  downstream  movement  of  the  tail 
of  the  vortex  patch;  see,  especially,  tin'  curves  marked  12  through  15. 

Figure  lb  shows  how  the  steady,  spanwise-periodic  vortex  pattern  gets  established  as  the 
hump  is  brought  to  its  final  shape  over  a  period  of  time.  Again,  one  may  note  tin*  monotonic 
behaviour  of  the  displacement  function  A(X.Z  =  0.  /  )  in  the  region  downstream  of  the  hump 
at  each  instant  ut  time.  In  the  vicinity  of  the  hump,  again,  tin*  displacement  function  closely 
resembles  the  instantaneous  hump  shape,  as  in  the  case  of  Fig.  la.  For  small  A.  this  latter 


behaviour  ran  also  be  predicted  using  the  closed  form  solution  (2.(i)  above,  which  further  shows 
that  for  A  <<  1,  the  wall-shear  perturbations  l  "'(0)  and  1F'(0)  an’  considerably  smaller  than 
the  displacement-thickness  perturbation,  being  of  0(.Y,/,!)  and  (9(A1//!),  respectively,  relative 
to  .1.  Flirt hermoie.  the  viscous  layer  begins  to  grow  at  the  usual  rate  of  Y  ~  A  1,1  near  A  <<  1. 
thus  suggest  ing  a  furt  her  split  with  in  t  he  lower  deck  corresponding  to  )  ~  A  A A  and  )  =0(1). 
The  wall-shear  scalings  for  small  A  show  that  the  flow  near  the  front  end  of  the  obstacle  is 
nearly  in  the  spanwise  direction:  however,  contrary  to  what  one’s  intuition  might  suggest,  this 
flow  is  found  to  be  converging  towards  the  obstacle  instead  of  spreading  out  along  the  spanwise 
direction.  A  physical  explanation  for  this  observation  is  provided  in  Section  3.2  below.  The 
unstable.  (Jbrtler  vortex  part  of  the  steady-state  solution  can  be  shown  to  be 

,1  =  A,  =  (-.W(0))3/s.  (3.2) 

■  > 

which  is  shown  via  symbols  in  Fig.  (  lb).  It  may  be  seen  1  ho.  the  numerical  solution  is  in 
excellent  agreement  with  the  instability-mode  part  of  it  for  nearly  all  locations  downstream  of 
the  hump.  Since’  tin-  algeoraically  decaying  part  of  the  solution  appears  to  be  dominated  so 
quickly  by  the-  unstable  eigenfunction  part,  this  suggests  that  experimental  measurements  of 
Cortler  vortices  need  not  only  be  carried  out  far  downstream  of  the  region  when’  the  vortices 
are  induced. 

Finally,  let  us  briefly  consider  the  inviscid  limit  corresponding  to  the  dispersion  relation 
( 2 . -S a ) .  It  is  easy  to  show  that  the  corresponding  differential  equation  for.  say,  the  displacement 
fi met  ion  .1  s  given  bv 

().x(c)x  +  \k-\(h)A-  k2A  =  Q.  (.3.3) 

which  is  hyperbolic  in  the  .V  -  T  space  on  the  basis  of  the  standard  classification  of  second 
order  partial  differential  equations.  The  two  families  of  characteristic  curves  corresponding  to 
equation  (3.3)  are  given  by 

7  =  constant  and  A  —  |F| /’ =  constant.  (3.1u.  M 

respectively.  This  suggests  that  the  Cauchy  problem  should  be  posed  at  some  initial  station 
.V  =  0,  and  if  the  Cauchy  data  has  a  finite  nonzero  limit  as  T  — ►  oc.  on<  duly  obtains  the  unstable 
Cbrtler  vortex  growing  ;n  the  downstieam  direction  as  part  of  the  causal  steady  solution.  In 
contrast,  if  one  traces  the  trajectories  of  both  roots  of  the  dispersion  relation  (2.Sa)  in  the 
complex  A  plane  as  rr  moves  from  positive  infinity  to  zero  along  the  real  axis  (fig.  o).  tin’ll 
application  of  the  Briggs  Bers  criterion  would  imply  that  the  root,  corresponding  to  the  unstable 
vortex  ought  to  appear  on  the  upstream  side.  The  reason  behind  this  failure  appears  to  be 
linked  to  the  infinite-speed  characlerist ic  (3.4a).  which  leads  to  a  solution  that  does  not  decay 
rapidlv  enough  in  A  at  anv  instant  of  time,  thereby  disallowing  the  usage  of  the  imaginary  A 
axis  as  the  inversion  contour  tor  a  two-sided  Laplace  transform  in  A  even  for  small  values  of  7  . 
Since  the  Briggs-  Bers  criterion  is  based  on  t  he  assn  m  pi  ion  of  a  spat  ially  compact  dist  nr  fiance  at 
small  enough  times,  it  seems  reasonable  to  expect  that  that  it  may  not  provide  reliable  results 
in  this  particular  case.  This  peculiar  behaviour  of  the  instability  in  question  is  also  reflected 
in  the  dispersion  relation  (2.Sa).  which  predicts  purely  dispersive  behaviour  for  real  streaimvise 
u  a  ven  u  m  bers  ( i  ,e.  neutral  behaviour  m  the  sense  of  a  t  emporal  inst  ability ).  but  predict  s  unstable 
roots  for  real  frequencies  (i.e.  spatial  instability),  in  particular,  the  stationary  (iortler  vorlex  at 

a.’  =  0. 


it) 


3.2  Unsteady  and  Steady  Disturbance  Patterns  due  to  Localized 
Surface  Irregularities 

In  view  of  the  order  of  magnitude  increase  in  the  computational  effort  required  for  this  rase, 
only  one  specific  situation,  corresponding  to 

F(A\  Z,T)  =  (1  -  co.s(2tt.V))  (1  -  co.s(4x7’))  c~^9  R(X,  0, 1 )  ft(7\  0, 0.5).  (3.5) 

was  considered  for  the  purpose  of  numerically  calculating  the  disturbance  pattern  due  to  an 
unsteady,  localized  wall  hump.  Here,  the  Gaussian  drop-off  in  the  spamvise  spectrum  was  chosen 
in  order  to  offset  the  ill-posedness  of  the  steady  dispersion  relationship  (2.9b).  The  contours  of 
constant  (instantaneous)  displacement  in  the  A'  —  Z  plane  for  this  case  have  been  shown  for  a  few 
selected  values  of  time  in  Figs.  (6a-e).  Due  to  the  spamvise  symmetry  of  the  surface  obstacle, 
only  positive  values  of  the  spamvise  coordinate  Z  have  been  included  in  the  plot.  Results  for 
negative  values  of  Z  can  be  obtained  through  a  reflection  of  the  contour  about  the  horizontal 
(A  )  axis.  Figures  6a-e  clearly  indicate  that  a  vortex  pattern  is  established  downstream  of  the 
source,  and  that  as  time  passes,  this  pattern  is  convected  downstream,  also  gaining  in  amplitude 
at  the  same  time.  This  further  demonstrates  the  convective  nature  of  Gortler  instability  in  the 
case  of  disturbances  with  a.  continuous  spectrum  of  spamvise  wavenumbers.  Thus,  it  would  be 
possible  to  study  the  disturbances  produced  by  a  steady  wall  hump  by  a  direct  solution  of  the 
steady  form  of  the  governing  equations. 

For  the  purpose  of  the  steady  calculation,  we  considered  the  hum])  shape  given  by 

F( A\  Z)  =  (  1  -  ro.s(2jr.Y))  c~z2/9  ft( .A,  0. 1 ).  (3.6) 

The  constant  displacement  contours  corresponding  to  the  disturbance  pattern  produced  by  this 
hum])  are  shown  in  Fig.  7a.  One  may  observe  that  the  ronforal  ellipses  near  the  inflow  boundary 
conform  rather  closely  with  the  obstacle  shape,  similar  to  that  in  the  spanwise-periodic  case 
discussed  in  Section  3.1  above.  Thus,  over  a  major  part  of  the  hump,  the  boundary  layer  is 
simply  lifted  up  by  the  amount  of  the  local  obstacle  height,  although  in  the  downward  sloping 
portion  of  the  hump,  there  are  regions  of  small  negative  displacement  at  sufficiently  large  values 
of  t  he  spamvise  coordinate  Z .  Figure  7a  also  indicates  the  presence  of  a  st ..ream wise  corridor  just 
downst  ream  of  the  hump  (corresponding  to  curve  6  closest  to  the  X-axis),  which  separates  the 
region  of  large  (positive)  displacement  from  the  relatively  less  disturbed  region  outside.  This 
feature  was  also  noted  by  Smith  ft  nl  (1977)  in  the  solution  to  a,  three-dimensional  triple  deck 
problem.  The  most  striking  feature  in  Fig.  7a,  however,  is  the  gradual  emergence  of  a  vortex 
structure  downstream  of  the  hump,  which  corresponds  to  alternating  regions  of  positive  and 
negat  ive  displacement,  separated  by  the  contours  marked  as  “6  \  At  each  streaniwi.se  stat  ion, 
the  innermost  vortex  is  the  strongest  one,  with  the  amplitude  falling  off  rather  rapidly  away 
from  the  centreline,  thus  leading  to  a  vortex  pattern  with  a  roughly  parabolic  shape'.  The 
orientations  of  the  centrelines  of  different  vortices,  corresponding  to  zero  displacement  contours 
(marked  as  “6  in  the  figure),  indicate  that  all  of  the  vortices  are  present  at  locations  just 
downstream  of  the  hump;  however,  since  their  amplitudes  are  inversely  proportional  to  the 
distance  away  from  the  axis  (Z  =  0),  the  more  distant  vortices  become  noticeably  large  only 
at  increasingly  larger  distances,  thereby  creating  the  illusion  that  additional  vortices  are  being 
created  as  one  moves  downstream.  The  qualitative  resemblance  between  the  vortex  patch  of 


Fig.  (7a)  ami  that  computed  numerically  for  the  case  of  (,'  =  0(1)  and  A-  =  0(1),  Denier.  Hall 
and  Seddougui  (1990).  where  no  asymptotic  approximations  were  involved,  indicates  that  the 
presence  of  spanwise  diffusion  terms  in  the  governing  equations  is  not  a  prerequisite  to  obtain  a 
slowly  spreading  pattern  of  st Teamwise  vortices. 

The  contours  of  constant  axial  and  spanwise  shear  stresses  (not  shown  here)  also  bear  some 
similarities  to  t  he  t  riple  deck  sit  nat  ion,  excepting  t  he  lack  of  upst  ream  influence,  and  t  he  presence 
of  the  dottier  vortex  pattern.  However,  an  examination  of  the  secondary  flow  within  the  X  —  Z 
plane  reveals  an  important  difference  between  the  disturbance  patterns  in  these  two  cases. 
Specifically,  the  results  of  Smith  tt  al  (1977)  show  that  on  the  forward  face  of  the  hump,  the 
fluid  is  pushed  out  in  the  spanwise  direction  within  the  lower  deck,  although  it  converges  again 
towards  t  he  obstacle  inside  t  he  middle  deck,  thus  creat  ing  a  recirculating  secondary-flow  pattern. 
This  pattern  is  reversed  on  the  backward  face  of  the  hump,  whereby  the  fluid  is  drawn  together 
in  the  lower  deck,  but  pushed  out  within  the  middle-deck.  Contrary  to  this,  as  well  as  to  one's 
intuitive  expectations,  it  is  found  that  in  the  case  of  an  obstacle  over  a  concave  surface,  the  fluid 
moves  towards  t  he  obstacle  on  bot  h  t  he  forward,  as  well  as  rear  fact's  of  the  hump.  This  spanwise 
convergence  of  the  fluid,  which  continues  across  the  middle  deck  as  well,  is  accompanied  by  an 
upwelling  near  the  axis  (Z  —  0).  thereby  setting  up  an  outflux  into  the  upper  deck,  where  the 
fluid  is  pushed  out  again:  see  Fig.  (7b).  The  counterrotating  vortices  generated  in  this  manner 
become  especially  prominent  farther  downstream  of  the  obstacle,  and  one  observes  the  distinct 
pattern  of  alternate  regions  of  upwelling  and  downwelling.  with  the  strength  of  this  pattern 
decreasing  quite  rapidly  away  from  the  axis. 

l  liese  differences  between  the  two  problems  can  be  easily  traced  to  the  nature  of  the  re¬ 
spective  relationships  between  the  pressure  and  the  displacement-thickness  perturbations.  In 
the  case  of  an  obstacle  with  triple-deck  scalings  over  a  flat  surface,  the  distribution  of  pressure 
is  controlled  by  the  streamwise  slope'  of  the'  displacement  function,  whereas  in  the  case  of  a 
streamwise  elongated  hump  over  a  concave  wall,  the  pressure  is  equal  to  the  negative  of  the 
displacement,  as  seen  from  the  viscous  limit  of  the  interactive  relationship  (2.. 'If).  Therefore,  the 
decreasing  displacement  away  from  the  centreline  in  Fig.  7a  translates  into  a  negative  pressure 
perturbation  of  decreasing  magnitude  away  from  the  axis.  Z  =  0.  In  other  words,  the  displace¬ 
ment  pattern  of  Fig.  7a  induces  a  favorable  pressure  gradient  towards  the  centreline,  which 
would  explain  the  sink-like  effect  the  obstacle  has  on  the  flow  within  the  boundary  layer. 

4  Conclusions 

We  have  investigated  the  unsteady  spatial  evolution  of  long  wavelength  (Jdrtler  vortices  in  a 
boundary  layer.  The  viscous-inviscid  interactive  regime  we  have  considered  is  appropriate  to 
situations  when  the  (Jdrtler  number  is  large  but  any  perturbations  in  the  flow  have  spanwise1 
wavelength  much  larger  than  the  boundary  layer  thickness.  In  particular  it  follows  that  our 
work  is  relevant  to  flows  near  the  leading  edge1  of  highly  curved  boundary  layers  of  the  type  one 
might  find  for  example  on  a  turbine  blade. 

We  have  shown  that  t  he  ( Jdrt  ler  mechanism  in  t  his  regime1  is  a  convective  rat  her  than  absolute 
instability.  This  result  contraelicts  previous  work  on  the1  (Jdrtler  problem  in  ibis  part  of  the 
(Jdrt lor- wavenumber  space.  Perhaps  the  main  conclusion  to  be  drawn  from  our  work  is  that 
in  general  it  is  sufficient  to  study  only  the  steady  (Jdrtler  problem.  However  we*  must  bear  in 


mind  that  our  analysis  is  restricted  to  the  small  wavenumber  regime  so  that  il  is  always  possible 
that  absolute  instabilities  might  occur  at  larger  wavenumbers.  The  fact,  that  experimental 
observations  all  indicate  that  Gortler  vortices  are  convectively  unstable,  suggests  that  the  latter 
scenario  is  unlikely. 
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Fig.  2  Spatial  root  of  the  steady  dispersion  relation  (2.7b) 
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Spatial  instability  properties  of  the  unsteady  dispersion  relation  (2.9a) 
for  k  =  1.0. 
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Fig.  4  Displacement  function  A(A',  Z,  T)  along  the  vortex  boundary  Z  =  0  as 
a  function  of  _Y  at  selected  values  of  time 

(a)  F(A,Z,T)  =  {1-cos(2kX))  (1-cos(4ttT))  cos(2nZ)  R(X,Q,l)R(T, 0,0.5) 
Curves  1  through  15  correspond  to  T=0.125,  0.250,  0.375,  0.5,  0.75, 

1.0,  1.25,  1.5,  1.75,  2.0,  2.25,  3.5,  6.0,  12.25,  and  18.5,  respectively. 


Fig.  4  Displacement  function  A( X,  Z,  T)  along  the  vortex  boundary  Z 
a  function  of  X  at  selected  values  of  time 


(1»)  F(X,  Z,  T)  =  (1  -  cos(2nX ))  (1  —  e 


T'/4))cos(2nZ)  Ii(X,  0, 1 ) 


Fig.  7a  Contours  of  constant  displacement  perturbation  for  the  steady  hump 


(3.6) 
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Fig.  7b  Sketch  of  the  secondary-flow  pattern  due  to  a  wall  hump  on  a  concave 
surface 
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